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SUMMARY 


A comparison is presented between the experimentally measured and 
theoretically calculated (by the method of NACA TN 3067 ) rates of roll 
of two rocket-propelled models with flexible rectangular wings. The 
comparisons show that although there are large aeroelastic losses in 
rolling rate, the theory predicts the actual rate of roll accurately. 


In reference 1, a method is presented for calculating the aero- 
elastic effects at supersonic speeds on the rolling behavior of air- 
craft with flexible rectangular wings . The method employs structural 
influence coefficients to determine the deformations of the wings and 
linearized supersonic lifting-surface theory to find the airloads . 

The purpose of the present paper is to assess the accuracy of the 
method of reference 1 by comparing its predictions with some experimental 
data presented in reference 2 for two, rocket -powered test models. Com- 
parisons between theory and experiment are given in the form of plots of 
rolling rate against Mach number. 


INTRODUCTION 


SYMBOLS 


D 


local flexural stiffness 



E 


Young's modulus of elasticity 


G 


shear modulus of elasticity 
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GJ 

r c /2 t 5( x ) 

elementary torsional stiffness, G / — - — dx 

J -c/2 5 

dGm(y,n) 

dy 

structural rate -of -twist influence function which results 
from a unit concentrated torque 

I 

moment of inertia of beam that represents flange effect 
of aileron 

K(x) 

stiffness of root springs 

L(y) 

aerodynamic load per unit span, positive upward 

M 

free -stream Mach number 

M(y) 

aerodynamic moment, per unit span, about midchord, 
positive in positive twist direction 

?h 

static pressure at altitude 

Po 

standard static pressure at sea level 

Q(y) 

aerodynamic moment, per unit span, about elastic axis, 
positive in positive twist direction 

V 

free -stream velocity 

w(y) 

deflection of midchord line of wing, positive upwards 

a 

ratio between fuselage radius and exposed wing semispan 

b 

total wing span, 2(aZ + l) 

c 

wing chord 

c a 

aileron chord 

e 

distance measured forward from midchord to elastic axis, 
expressed as fraction of chord 

l 

exposed wing semispan 

P 

rolling velocity 

pb 
2 V 

tangent of wing -tip helix angle, positive for counter- 
clockwise roll when viewed from behind 
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s 

ratio of root flexibility to wing flexibility 

t 

thickness of wing cross section 

w(x,y) 

local deflection of wing, positive upwards 


coordinate system 

X 

distance from midchord to location of beam that represents 
flange effect of aileron 

e(y) 

angle of twist of wing, leading edge up 

9 r 

twist at mid exposed span which results from a unit 
torque at the tip 

6 

aileron deflection, positive down 

A 

nondimens ional parameter involved in the rate-of -twist 
influence function (see eq. (A10)) 


Poisson's ratio 

9 

(pb/2V)p 

rolling effectiveness, 

(pb/2V) R 


Subscripts : 


F 

flexible wing 

R 

rigid wing 

rev 

aileron reversal 


DESCRIPTION OF THE MODELS 


The two models investigated (the last two in table I of ref. 2) were 
essentially the same in size and shape, the only important difference 
being that one had aluminum wings and the other had steel. Both models 
had three rectangular wings equally spaced around the rear portion of a 
long cylindrical body. (See fig. 1.) The wings were uniform in the 
spanwise direction and had NACA 65AOO3 airfoil sections. The full-span 
trailing-edge ailerons were formed by bending the wing along the 
80-percent-chord line. 
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The pertinent dimensions of the two models are given in table I . 
Included in table I are the experimentally determined values of 0 r , 

the twist at the midspan due to a unit torque at the tip. These values, 
which were obtained from reference 2, are used as an aid in determining 
the structural characteristics of the wings. The ordinates for the NACA 
65AOO3 airfoil, obtained from reference 3, are also included in this 
table . 


THEORETICAL RESULTS 


The experimental data in reference 2 are given in the form of plots 
of rolling rate pb/2V against Mach number M. Also given are the vari- 
ations with Mach number of altitude (specified in the form of static 
pressure) during the flights. It is desired to calculate theoretically 
the variation of pb/2V with Mach number for the prescribed altitude 
variation for each of the two models. In order to do this, the compu- 
tational procedure outlined in reference 1 is followed. 

In reference 1, several alternative calculation schemes were 
described. The particular one used herein is the same as that used in 
the analysis of the example configuration in reference 1. This approach 
is exemplified by the matrix equations (29) or (31) of that report. 

These equations are written in terms of the rate of twist d0/dy rather 
than the twist 0 itself and make use of an interpolation procedure to 
reduce the number of degrees of freedom involved. For the present prob- 
lem the structural ingredients of these matrix equations - the rate-of- 
twist influence functions - are derived in the appendix of this report 
by the application of an approximate plate theory and are tabulated for 
both models in table II j the aerodynamic ingredients - the various 
indicial loads and the loads due to roll and aileron deflection - are 
obtained from reference 1. (The assumption is made that the loads on 
each wing of the three -winged aircraft considered in this paper are the 
same as those resulting from the two-winged configuration considered in 
ref. 1.) It should be noted that the analysis in the appendix indicates 
the existence of an "elastic axis," a line along which loads can be 
placed without producing any appreciable twist. Accordingly, equa- 
tions (29) and (31) of reference 1 are modified as suggested therein to 
take advantage of this elastic axis . 

By using these modified equations, then, and by following the sug- 
gested computational procedure, the theoretical results presented herein 
were obtained. The results for aileron reversal appear in terms 
of /P Ip ) , the ratio between the static pressure at which the ailerons 

l hi o) rey > 

reverse and the standard sea-level static pressure . In figure 2 is shown 
the variation with Mach number of this ratio for both the aluminum and 
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steel wings . 


A calculation of the rolling effectiveness 



for other values of static pressure at several different Mach numbers 
showed that the variation of cp with P^Pq was almost exactly linear. 


Therefore , the values of cp for any pressure ratio at a given Mach num- 
ber can be deduced from the value of pP^y/P^ for that Mach number . 


Thus, for either model 



( 1 ) 


The rolling rate pb/2V 
multiplying cp by (pb/2V) R . 


for the flexible wing can be obtained by 
Consequently, 




( 2 ) 


The variation of 



with Mach number has been found by the 


method of reference 1; these values were used in conjunction with the 
information in figure 2 and the plots of the actual flight values 
of P h /P o against Mach number from reference 2 to obtain the theoretical 

values of ^/5 shown in figure 3 for the two models. Also shown in 
figure 3 is the theoretical for the rigid wing. It should be 

remarked that measured values of P h /P Q were lacking for Mach numbers 

higher than 1.4 for the steel wing and 1.8 for the aluminum wing; accord- 
ingly, the curves for the flexible wings have been stopped at these values. 


COMPARISONS AND DISCUSSION 


The variation of experimental |^/s with Mach number, as obtained 

from reference 2, is shown in figure 3 f° r comparison. The estimated 
experimental rigid rate of roll, obtained by extrapolation from the 
flexible data by assuming a linear variation of pb/2V with the param- 
eter P R 0 r , is also shown in this figure. 
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From figure 3 it can be seen that although there is a large loss 
of rolling effectiveness due to aeroelasticity, the theory does a good 
job of predicting the actual flexible rate of roll. 


CONCLUSION 


A comparison between theoretically and experimentally determined 
rates of roll for two rocket-propelled models with flexible rectangular 
shows that the method of NACA TN 30^1 rs capable of yielding 
accurate predictions for the aeroelastic effects on the roll of super- 
sonic aircraft with rectangular wings. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 11, 195^* 
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APPENDIX 

DERIVATION OF TORSIONAL INFLUENCE COEFFICIENTS 


In order to solve the aeroelastic rolling problem it is necessary 
to compute the torsional influence coefficients for the wing. These 
quantities are derived in this appendix by means of the same approach 
as that of reference 4; that is, the deflections are assumed to be lin- 
ear in the chordwise direction and the principle of mi niminn potential 
energy is employed. 

The structure under consideration is shown in figure 4(a). It con- 
sists of a solid plate that is uniform in the spanwise direction with a 
bent-up aileron and a more -or -less -rigid attachment to the model body. 

Two factors prevent the analysis of this structure directly by the method 
of reference 4: (l) The bend along the aileron hinge line produces a 

flange effect so that flat-plate theory cannot be used; ( 2 ) the root of 
the wing cannot be considered to be perfectly clamped. The analysis is 
therefore performed for the equivalent structure shown in figure 4(b). 

In this figure, the aileron has been unbent and the flange effect has 
been represented by a beam. The moment of inertia of the beam is assumed 
to be equal to the difference between the moments of inertia of the bent 
and unbent aileron and the beam is located at the centroid of the 
difference -in-moment-of -inertia distribution. Thus, 



Also, in order to represent the effect of incomplete root clamping, the 
plate-beam combination is assumed to be mounted on springs which prevent 
displacement of the root but permit non-zero slopes in the spanwise 
direction. In the analysis to follow, the shape of the stiffness 
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distribution of these springs, which, for the present, is indeterminate, 
is assumed to be such that the resulting equations exhibit their sim- 
plest possible form; the absolute magnitude of the spring stiffness for 
each of the two models is then selected so that the theoretical twist at 
the midspan due to a unit torque at the tip matches the experimental value 
in table I. 


The potential energy of this equivalent structure subjected to the 
distributed lateral load p(x,y) is 




pi pc/ 2 

0 -c/2 


D(x) 


!$)“•( 


•« - & 


>dxdy + 


El 


r l 

8 2 w(x,y) 

>0 

dy 2 


y 


i r /2 

2 ^-c/2 


K(x) 


(x,0) 

oy 


>* dx - 


'0 


>1 p c/2 

-c/2 


dy + 

p(x,y)w(x,y) dxdy 

(A3) 


where w(x,y) is the deflection, D(x) = 


Et' 


2 ) 


is the plate stiff - 


12(1 - p 

ness, and K(x) is the stiffness distribution of the root springs. 


In accordance with the assumption of linear chordwise deformations. 


let 


w(x,y) = W(y ) - x0 (y ) 


(A4) 


Then, upon performing the integration with respect to x. 
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+ 



El) (W") 2 - 

(e ") 2 + 2(1 

- 2k 2 W'0' + 


2 (a.2 + EHr)w"0" + 


m2 


- M- ) a 1 (0 ' ) 

^3(0 ' )^1 

-*y=o 


dy + 



|L(y)W(y) + M(y) 0 (yj]dy 


where 


and 



(A5) 


are, respectively, the section lift and moment about the midchord. 

Minimization of the potential energy yields the following differ- 
ential. equations, 
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(a ± + El)w IV - (a 2 + EIx)e IV = L 
- (a 2 + EIx)w IV + (a 5 + EIx 2 ) 0 IV - 2(1 - ia)a 1 9” 


> 


= M 




(a6) 


and boundary conditions. 


w( 0 ) = 9(0) = 0 

(a. ± + EI)W"(0) - (ag + E33c)9"(0) - *^'(0) + k 2 9’(0) = 0 

(a 2 + Elx)w"(0) - (a 5 + EIx 2 ^"^) - k 2 W'(0) + k 5 9'(0) = 0 
(ai + El)w"(z) - (a 2 + EIx)9"(z) = 0 

(a. 2 + EIx)w"(z) - (a^ + EIx^9"(z) = 0 
+ El) W" * ( Z ) - (a 2 + EIx)q " ' ( Z ) = 0 

(a 2 + EIx W"'(Z) - (a^ + EIx 2 9"'(Z) + 2(l - ti)a 1 9 , (2) = 0 

For purposes of aeroelastic calculations, only the twist 9 is 
important; therefore, it is desirable to eliminate W from equations (A6) 
and the accompanying boundary conditions. In order to do this easily, 
it is convenient to take advantage of the freedom of choice of the shape 
of the root spring stiffness distribution and assume that 


k^ k 2 k^ 

a-L + El a 2 + EIx a + El5 j2 

Since, as has been mentioned, the overall spring stiffness is to be 
selected by duplicating experimental values of twist, this l im itation 
on the distribution should have negligible influence on the desired 
aeroelastic results . 
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After eliminating W, the resulting equations are 


_2 (ao + EIx) 2 

az + EIx - !— 

■> a]_ + El 


e 17 - 2(1 - n^e" = m + 


ap + EIx . 

-2 — L A7 

a^ + El 


0 ( 0 ) = 0 


9 "( 0 ) = 


aj_ + El 


0'(O) 


0"(z) = o 


+ BETS - 


a^ + El 


0 1 ' ' ( 2 ) - 2(1 - M)a 1 0'(Z) = 0 


Integrating equation (A7) once and rewriting yields 


0'" - % 0 1 =-4?T7 \ l Q(n) dn 

i l * ^y 


(A8) 


where the remaining boundary conditions are 

">v 

0 ( 0 ) = 0 


0"(O) = -i- 0'(o) 

sZ 


0 "( Z ) = 0 


(A9) 
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and where 


A = 


2(1 - ix)^ 


—\2 


a, + EIx 2 - 
5 


( a -2 + Elxj 
aj_ + El 


(A10) 


GJ = 2(1 - d)a 1 
Q(y) = M(y) - ecL(y) 

&2 + EIx 
c (a.^ + El) 

a-i + El 

s = — 


The parameter GJ can be written 


GJ = G 


'° /2 tS(x) 


-c/2 


dx 


(All) 

(A12) 

(A13) 


(Aik) 


This quantity then is merely the elementary torsional sitffness for 
a thin cross section. The quantity Q is the total moment about an 
"elastic axis" located e chords ahead of the midchord. The term 
"elastic axis" is justified in this case because, as can be seen from 
equation (A8) , loads applied at x = -ec anywhere along the span would 
produce no twist. The parameter s expresses the ratio between the 
overall stiffnesses of the wing and of the root springs. 


In view of the existence of an elastic axis for this configuration, 
only the influence functions due to torque need be obtained. More spe- 

dfi 

cifically, it is desired to calculate the rate of twist — at sta- 

dy 

tion y caused by the application of a unit torque at station q. This 


quantity, designated 
with the appropriate 


^2S.(y,T]), can be obtained by solving equation 

dy 

substitution for Q(y) and is given by 


(A8) 
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% y ,*) . L c ° sh 4 ~ I) 

6 y GJ cosh A + sA sinh A 


"A 


1 + cosh A — + 

l 


sA sinh A ^ 


dy 


(y,n) = -± 


GJ cosh A + sA sinh A 


cosh A - cosh A 1 - *- - ^ 


(y > n) 


sinh A (1 - jj sinh A ~ + 


sA 


sinh A - sinh A (1 - ^ j cosh A y 


(y < tj) 

-/ 

(A15) 


All the quantities in these equations except s can be found 
directly from the geometry and material properties of the wing. The 
parameter s can be evaluated by equating the theoretical and experi- 
mental values of the twist at midspan due to a torque at the tip. Such 
a process yields the following formula: 



coth A 
A 


(Al6) 


It should be noted that by determining s in this manner, the 
torsional characteristics of the wing are duplicated closely. If the 
value of s were taken to be zero, the root would be completely clamped 
and the wing would be too stiff. If, on the other hand, a value of s 
equal to infinity were used, the root would be completely free to warp 
and the wing would be too flexible. For the aluminum wing, for example, 
taking s to be zero and infinity would yield, respectively, 0 r equal 

O .956 x 10 _1+ and 1.262 x 10 _lf radians in. -lb, whereas the actual experi- 
mental value is 1.175 radians in. -lb. From these values, it can be seen 
that appreciable error could result from an improper selection of s. 

In order to evaluate the rate-of -twist influence functions, the 
quantities I and x (eqs . (Al) and (A2)) were computed for the nominal 
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aileron deflection of 5°; in this calculation, the thickness t(x) was 
assumed to vary linearly from its value at the hinge line (see table I) 
to zero at the trailing edge. The values of a 1? a 2 , and a^ (defined 

following eq. (A5)) were computed from the airfoil profile data in 
table I; in this computation, the necessary integrations were performed 
numerically by using Simpson's rule. The quantities 'A, GJ, and e 
were then evaluated. (See eqs . (A10), (All), and (A13)») Finally, equa 
tion (Al6) was used in conjunction with the experimentally determined 
values of 0 r to find s. A tabulation of the pertinent parameters for 

each of the two models follows : 


Steel Aluminum 

0.0479 0.0485 


GJ, lb-in. 2 117,500 42,100 

A 8.4879 8.3138 

O.O891 0.289 


The values of dG m jd 

computed for 0 < — < 1 

~L 


y for the aluminum and steel wings have been 

in steps of 0.2 and 0 < — < 1 in steps of 0.1 

l 


The results are given in table II in matrix form. These matrixes are 
used directly for the computation of the aeroelastic effects on roll by 
the method of reference 1. 
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TABLE I.- DESCRIPTION OF MODELS 


(a) Basic information 


Parameter 

Steel model 

Aluminum model 

l, in 

10.61 

10.61 

c , in 

7.O7 

7.07 

a 

0.236 

0.236 

c a/ c 

0.2 

0.2 

E, psi 

29 X 10 6 

10.6 X 10 6 

M- 

0.300 

0.333 

0 r , radian/in . -lb . . 

0.392 X 10 _ 4 

a 1.175 x 10 -4 

5, deg 

5-5 

4.9 


a Two values of 0 r for the aluminum model were given in 
reference 2. The value herein is the correct one. 


(h) NACA 65AOO3 airfoil ordinates 


x/c 

t/c 

x/c 

t/c 

x/c 

t/c 

x/c 

t/c 

-0.500 

-.495 

-.4925 

-.4875 

-.475 

-.45 

-.425 

-.400 

0 

.00464 

.00564 

.00718 

.00982 

.01314 

.01592 

.01824 

-0.350 

-.300 

-.250 

-.200 

-.150 

-.100 

-.050 

0.0 

0.02194 

.02474 

.02688 

.02842 

.02946 

.02996 

.02992 

.02926 

0.050 

.100 

.150 

.200 

.250 

.300 

.350 

.400 

0.02794 

.02606 

.02364 

.02088 

.OI776 

.01438 

.01090 

.OO728 

0.450 

.500 

0.00370 

.00014 
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TABLE II.- RATE-OF -TWIST INFLUENCE-FUNCTION MATRIXES 


(a) Steel Wing, GJ = 117,500 lb-ln. 2 




0 

0.246330 

0.351742 

0.396851 

0.416156 

0.424417 

0.427954 

0.429471 

0.430127 

0.430426 

0.430595 



0 

0.115132 

0.398052 

0.682758 

0.804593 

0.856734 

0.879056 

0.888629 

0.892773 

0.894659 

0.895727 

sg m 

= i_ 

0 

0.021085 

0.072897 

0.195062 

0.481008 

0.767011 

0.889456 

0.941967 

0.964701 

0.975043 

0.980904 

Sy 

" G J 

0 

0.003865 

0.013364 

0.035760 

0.088182 

0.210716 

0.497071 

0.783847 

0.908006 

0.964492 

0.996499 



0 

0.000731 

0.002526 

0.006760 

0.016671 

0.039837 

0.093973 

0.220480 

0.516103 

0.918371 

0.999338 



o_ 

0.000259 

0.000895 

0.002396 

0.005908 

0.014117 

0.033301 

0.078132 

0.182891 

0.427698 

0.999765 


(b) Aluminum Wing, GJ = 42,100 lb-in. 2 




0 

0.398616 

0.572192 

0.647775 

0.680688 

0.695021 

0.701265 

0.703991 

0.705193 

0.705750 

0.706073 



0 

0.144975 

0.436856 

0.723321 

0.848060 

0.902382 

0.926048 

0.936378 

0.940933 

0.94307 

0.944267 


= _l_ 

0 

0.027491 

0.082838 

0.206553 

0.489192 

0.771643 

0.894695 

0.948408 

0.972092 

0.983085 

0.989432 

Sy 

~ GJ 

0 

0.005219 

0.015727 

0.039215 

0.092874 

0.215979 

0.498635 

0.781750 

0.906593 

0.964537 

0.997994 



0 

0.001024 

0.003085 

0.007693 

0.018220 

0.042370 

0.097820 

0.225155 

0.517578 

0.823168 

0.999606 



0 

0.000375 

0.001129 

O.OG2816 

0.006670 

0.015510 

0.035909 

0.082422 

0.189469 

0.435303 

0.999856 
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Figure 2.- Variation of pressure ratio at reversal with Mach number for 

the two models. 
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1.0 1.2 1.4 1,6 1.8 


Mach number, M 

Figure 3-- Variation with Mach number of the theoretical and experimental 

•rate of roll. 
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(a) Wing with bent-up aileron and flexible root attachment. 



A 



(d) Equivalent plate-beam combination mounted on root springs. 


Figure k .- 


Wing with bent-up aileron and equivalent plate-beam configuration. 
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